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ABSTRACT 



' We report the discovery by the WASP transit survey of a giant planet in a close orbit (0.0295 ± 0.0009 AU) around a moderately 

bright (V = 11.6, K = 10) G9 dwarf (0.89 ± 0.08 M Q , 0.84 ± 0.03 R e ) in the Southern constellation Eridanus. Thanks to high- 
precision follow-up photometry and spectroscopy obtained by the telescopes TRAPPIST and Euler, the mass and size of this planet, 
WASP-50b, are well constrained to 1.47 ± 0.09 M Jup and 1.15 ± 0.05 R ]up , respectively. The transit ephemeris is 2455558.6120 
(±0.0002) + N x 1 .955096 (±0.000005) HJD UTC • The size of the planet is consistent with basic models of irradiated giant planets. The 
chromospheric activity (log R' HK = -4.67) and rotational period (P l0t = 16.3 ± 0.5 days) of the host star suggest an age of 0.8 ± 0.4 Gy 
that is discrepant with a stellar-evolution estimate based on the measured stellar parameters (p, = 1 .48 ± 0. 10 p , T c g= 5400 ± 100 K, 
[Fe/H]=-0.12 ± 0.08) which favours an age of 7 ± 3.5 Gy. This discrepancy could be explained by the tidal and magnetic influence of 
the planet on the star, in good agreement with the observations that stars hosting hot Jupiters tend to show faster rotation and magnetic 
activity (Pont 2009; Hartman 2010). We measure a stellar inclination of 84+*, deg, disfavouring a high stellar obliquity. Thanks to its 
large irradiation and the relatively small size of its host star, WASP-50 b is a good target for occultation spectrophotometry, making it 
able to constrain the relationship between hot Jupiters' atmospheric thermal profiles and the chromospheric activity of their host stars 
proposed by Knutson et al. (2010). 



Key words, stars: planetary systems - star: individual: WASP-50 - techniques: photometric - techniques: radial velocities - techniques: 
spectroscopic 
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• 1. Introduction Jupiters also make it possible to study the physical response 

__ ' _ ., , , ... ... . of a giant planet to an irradiation orders of magnitude larger 

^ ■ While the Kepler space mission is pursuing its pioneering than fa solar m lanets (Burrows et al 2008; Fort et al 

, search for habitable terrestrial planets transiting solar-type stars 2010)> and ^ tQ ±e g grav i tationa i and magn etic fields 
m > ■ (Borucki et al. 201 1) ground-based, wide-field surveys continue &Q dose tQ ^ hofjt starg (CojTeia & Lagkar 2010; ch gt 
^ ; to detect short-period (i.e. a few days, or even less) transiting gi- aJ 201Q) Furthermore? these lanets could be able t0 influe nce 
KJS . ant planets at an increasing rate. These transiting hot Jupiters ±e ert ies of their t stars> notabl b modifying their 
^ ■ are important objects for the nascent field of comparative exo- ^ momentum budget and by inducing chromospheric ac- 
planetology. As they could not form in situ at such short dis- dy{ (Lanza 201Q) The atm neres of hot j iters can also 
tances from their parent stars, their orbital configurations pose fee smdied th hl Their generally large atmospheric scale 
an interesting challenge for theories of formation and evolu- hd hfa and ^ fr transks maximize the signa i_ to _ no i se 
tion of planetary systems (D Angelo et al. 2010; Lubow & Ida ratios (§NR) achievable on short timesC ales in transit absorp- 
2010). Their migrational history can be further constrained by tion spectrophotometry, while their large irradiation makes it 
the measurement of the orientation of their orbital axis relative iWe tQ ^ themal emission ffle of ^ d 
to the stellar spin axis (Winn 2010a), the measurements obtained ^ ^ h occultation . dept h measurements at different wave- 
so far suggesting that past dynamical interactions with a third , ^ (S 2010) Such atmo heric measurements have 
body could be driving the ir migration (Triaud et al. 2010). Hot made possible ±e smdy of ^ atmospheric composition of sev- 

~ , „ . I , . .„ „ . , eral planets (Seager & Deming 2010). They have also revealed 

Send offprint requests to: michael.gillon@ulg.ac.be , r , . , ,. 

* „, , . , . , . , . .... that most hot Jupiters harbour a high-altitude thermal inversion 

The photometric time-series used in this work are only available , , r , ■ ■ •,• 

in electronic form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr on thelr da y slde ' the Presence or absence of this inversion being 

(130.79.128.5) or via http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/ possibly correlated with the activity of the host star (Knutson 
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et al. 2010). Furthermore, occultation measurements can help in 
measuring their orbital eccentricity, an important parameter to 
assess their energy budget and tidal history (Gillon et al. 2010a). 
These planets are thus real 'exoplanetology laboratories', espe- 
cially if they orbit stars bright enough to make possible high- 
SNR follow-up studies. 

In that context, we announce here the discovery of a hot 
Jupiter, WASP-50b, that transits a moderately active, V=11.6, 
G9V star in the Southern constellation Eridanus. Section 2 
presents the discovery photometry obtained by the WASP tran- 
sit survey, and high-precision follow-up observations obtained 
from La SiUa ESO Observatory by the TRAPPIST and Euler 
telescopes to confirm the transits and planetary nature of WASP- 
50b, and to determine precisely its parameters. In Sect. 3, we 
present the spectroscopic determination of the stellar properties 
and the derivation of the system parameters through a combined 
analysis of the follow-up photometric and spectroscopic time- 
series. We discuss our results in Sect. 4. 



2. Observations 

2.1. WASP photometry 

The host star WASP-50 (1SWASPJ025445. 13-105353.0 = 
2MASS02544513-1053530 = USNO-B1.0 0791-0028223 = 
Tycho-2 5290-00462-1 = GSC1 5290-00462 = GSC2.3 
S2L3000257; 7=11.6, K=10.0) was observed by the Northern 
and Southern telescopes of the WASP survey (Pollacco et al. 
2006) during the 2008 and 2009 observing seasons, covering the 
intervals 2008 Sep 16 to 2009 January 3 and 2009 August 30 to 
2009 December 19. The 19731 pipeline-processed photometric 
measurements were de-trended and searched for transits using 
the methods described by Collier Cameron et al. (2006). The se- 
lection process (Collier Cameron et al. 2007) found WASP-50 
as a high priority candidate showing a periodic transit-like sig- 
nature with a period of 1.955 days. A total of 26 transits are 
observed in the data. Fig. 1 (top panel) presents the WASP pho- 
tometry folded on the best-fitting transit ephemeris. 

We analyzed the WASP lightcurve to determine whether it 
shows periodic modulation due to the combination of magnetic 
activity and the rotation of the star. The value of v sin i m derived 
in section 3.1 together with the estimated stellar radius implies 
a rotation period of 16-24 days, assuming that the rotation axis 
of the star is approximately aligned with the orbital axis of the 
planet. We subtracted the transit signal from the data and then 
used the sine-wave fitting method described in Maxted et al. 
(2011) to calculate the periodograms shown in Fig. 2. These 
are calculated over 4096 uniformly spaced frequencies from to 
1.5 cycles/day. The false alarm probability (FAP) levels shown 
in these figures are calculated using a bootstrap Monte-Carlo 
method also described in Maxted et al. (2011). Variability due 
to star spots is not expected to be coherent on long timescales, 
as a consequence of the finite lifetime of star-spots and differen- 
tial rotation in the photosphere, so we analysed the two seasons 
of data for WASP-50 separately. 

Both seasons of data shows a significant periodic modulation 
with an amplitude of 4-5 milli-magnitudes, although at slightly 
different periods — 16.06 + 0.09 days and 16.65 ± 0.13 days. 
The FAP is smaller than 0.05% for season 1, and ~ 0.1% for 
season 2. We used a bootstrap Monte-Carlo method to estimate 
the error on the periods. We also calculated periodograms of two 
stars with similar magnitudes and colours to WASP-50 observed 
simultaneously with the same camera. These stars show signifi- 
cant periodic modulation at periods of 27 and 160 days and with 




Fig. 1. WASP photometry of WASP-50 folded on the best-fitting 
transit ephemeris from the transit search algorithm presented in 
Collier Cameron et al. (2006), and binned per 5 measurements. 



amplitudes of about 1 1 milli-magnitudes in the first season of 
data. One of the stars shows significant variability in the second 
season of data near a period of 1 day, the other shows no sig- 
nificant variability in this frequency range. Our interpretation of 
these results is that WASP-50 does show a periodic modulation 
of about 16.3 days, but that instrumental noise in the WASP data 
adds a systematic error of about ±0.5 days to this estimate of the 
rotation period. Combining the periods from both seasons and 
taking into account this systematic error, we conclude that the 
rotation period of WASP-50 is 16.3 + 0.5 days. One could argue 
that the distribution of star spots on the surface of a magneti- 
cally active star can produce periodic photometric variations at 
half of the true rotation period. This is not very likely in the case 
of WASP-50 because a star with a rotation period of 32 days is 
likely to have very low levels of magnetic activity, and the im- 
plied value of vsin/» (< 1.3 km s^ 1 ) would be inconsistent with 
the observed value given in Table 2. 

2.2. High-SNR transit photometry 

2.2.1 . TRAPPIST / + z-band photometry 

In order to confirm the origin of the transit signal on WASP- 
50 and to constrain thoroughly its parameters, three transits 
were observed at high-SNR with the new telescope TRAPPIST 
located at ESO La Silla Observatory in the Atacama Desert, 
Chile. TRAPPISlQ (77Mnsiting Planets and Planetes/mals 
Small Telescope; Gillon et al. 2011) is a 60cm robotic tele- 
scope dedicated to high-precision photometry of exoplanet host 
stars and to the photometric monitoring of bright comets. It was 
installed at La Silla in April 2010. Its commissioning ended 
in November 2010, and it has been operational since then. Its 
thermoelectrically-cooled camera is equipped with a 2K x 2K 
Fairchild 3041 CCD with a 22' x22' field of view (scale = 
0.65"/pixel). We observed three consecutive transits of WASP- 
50 on the nights of 2010 December 23, 25 and 27. We used the 
2x2 MHz read-out mode with lxl binning, resulting in a 
typical readout + overhead time and readout noise of 6.1 s and 
13.5 e~, respectively. A slight defocus was applied to the tele- 
scope to optimize the observation efficiency and to minimize 
pixel to pixel effects. We observed through a special 'I+z' fil- 
ter that has a transmittance >90% from 750 nm to beyond 1 100 
nm. The positions of the stars on the chip were maintained to 
within a few pixels over the course of the three runs, thanks to the 
'software guiding' system that regularly derives an astrometric 



see http://arachnos.astro.ulg.ac.be/Sci/Trappist 
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Fig. 2. Upper panels: periodograms for the WASP data from two seasons for WASP-50. The date range (JD-245000) is given in 
the title of each panel. Horizontal lines indicate false alarm probability levels FAP=0. 1,0.01,0.001. Lower panels: data folded and 
binned in 0.005 phase bins for the periods noted together with the date range in the plot title. 



solution to the most recently acquired image and sends pointing 
corrections to the mount if needed. 

Weather conditions were good during the three runs. The in- 
tegration time was 25 s for the first run and 20s for the other 
runs. After a standard pre-reduction (bias, dark, fiatfield correc- 
tion), the stellar fluxes were extracted from the images using the 
IRAF/DAOPHOtJ] aperture photometry software (Stetson, 1987). 
Several sets of reduction parameters were tested, and we kept 
the one giving the least scatter for stars of similar brightness to 
WASP-50. After a careful selection of reference stars, differen- 
tial photometry was then obtained. Fig. 3 shows the TRAPPIST 
transit light curves divided by the best-fitting baseline model, 
with the best-fitting transit model superimposed. 

2.2.2. Euler r-band photometry 

The transit of 25 Dec 2010 was also observed in the r Gunn filter 
with the EulerCam CCD camera at the 1 .2-m Euler Telescope at 
La Silla Observatory. This nitrogen-cooled camera has a 4k x 4k 
E2V CCD with a 15' x 15' field of view (scale=0.23"/pixel). 
The exposure time was 40 s, the readout + overhead time be- 
ing 25 s (readout noise ~4.5e~). Here too, a defocus was ap- 
plied to the telescope to optimize the observation efficiency and 
minimize pixel to pixel effects, while flat-field effects were fur- 
ther reduced by keeping the stars on the same pixels, thanks to 
a 'software guiding' system similar to TRAPPIST's. However, 
a residual drift of the target on the detector was observed for 
airmasses larger than 1.5. This positional drift did not result in 



2 IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 



any clear drift in the photometry itself. Furthermore, any po- 
tential systematic effect on the photometric baseline is included 
in our analysis under the form of a second-order time polyno- 
mial (see Sect. 3.2). The reduction was similar to the one of the 
TRAPPIST data, i.e. aperture photometry was extracted for all 
stars in the pre-reduced images, then differential photometry was 
obtained. The resulting light curve is shown on Fig. 3 (bottom). 

2.3. Spectroscopy and radial velocities 

Once WASP-50 was identified as a high priority candidate we 
gathered spectroscopic measurements with the C0RALIE spec- 
trograph mounted on Euler to confirm the planetary nature of 
the transiting body and obtain a mass measurement. 15 spectra 
were obtained from 2010 December 5 to 201 1 January 5 with an 
exposure time of 30 minutes. Radial velocities (RV) were com- 
puted by weighted cross-correlation (Baranne et al. 1996) with a 
numerical G2-spectral template (Table 1). 

The RV time-series shows variations that are consistent with 
a planetary-mass companion. A preliminary orbital analysis of 
the RVs resulted in a period and phase in excellent agreement 
with those deduced from the WASP transit detections (Fig. 4). 
Assuming a stellar mass M* = 0.91 ± 0.07 Mq (Sect. 3), the 
fitted semi-amplitude K = 260 + 5 ms _1 translates into a plan- 
etary mass slightly larger than Jupiter's, M p - 1.52 + 0.08 
Mj up . The resulting orbital eccentricity is consistent with zero, 
e = O.OI^qJJj. The scatter of the residuals is 10.8 ms _I , very 
close to the average RV error, 10.5 ms~'. 

To exclude the possibility that the RV signal originates from 
spectral line distortions caused by a blended eclipsing binary, 
we analyzed the C0RALIE cross-correlation functions using the 
line-bisector technique described in Queloz et al (2001). The bi- 
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Fig. 3. Left, from top to bottom: period-folded TRAPPIST transit light curves of WASP-50 for 2010 December 23, 25, and 27, and 
Euler transit light curve. These light curves are divided by the best-fitting baseline models (see details in Sec. 4), and the best-fitting 
transit models are superimposed in red (TRAPPIST) and green (Euler). Three light curves and models are shifted in the y-axis for 
the sake of clarity. Right: the residuals to the fit for each light curve. 



sector spans seem quite stable, their standard deviation being 
nearly equal to their average error (22 vs 21 ms~'). No evidence 
for a correlation between the RVs and the bisector spans was 
found (Fig. 5), the deduced correlation coefficient being -0.24. 
Considering not only this absence of correlation, but also the 
excellent agreement between the stellar densities inferred from 
the high-SNR transit photometry (Sect. 3.2) and from the spec- 
troscopic analysis (Sect. 3.1), we concluded that a giant planet 
transits WASP-50 every -1.96 days. 



3. Analysis 

3.1. Spectroscopic analysis - stellar properties 

The CORALIE spectra of WASP-50 were co-added to produce 
a single spectrum with a typical SNR of 100. The standard 
pipeline reduction products were used in the analysis. 

The analysis was performed using the methods given by 
Gillon et al. (2009a). The H a line was used to determine the 
effective temperature (r e jf), while the NaiD and Mgib lines 
were used as surface gravity (logg) diagnostics. The parameters 
obtained from the analysis are listed in Table 2. The elemen- 
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Fig. 4. Top: Euler/CORALIE RVs with the best-fitting Keplerian 
model superimposed in red. Bottom: residuals of the fit. 




RV [kms-1] 

Fig. 5. Bisector versus RV measured from the CORALIE spectra. 
We adopt uncertainties of twice the RV uncertainty for all bi- 
sector measurements. There is no correlation between these two 
parameters indicating the RV variations are not caused by stellar 
activity or line-of-sight binarity. 



HJD-2,450,000 


RV 


o~rv 


BS 




(km s" 1 ) 


(m s- 1 ) 


(km s" 1 ) 


5535.647447 


25.2667 


8.5 


-0.0341 


5536.647197 


25.7970 


8.5 


-0.0363 


5537.656077 


25.2840 


10.1 


0.0054 


5538.646866 


25.7690 


9.3 


-0.0554 


5541.573212 


25.2571 


8.8 


0.0006 


5541.735749 


25.3342 


11.0 


-0.0410 


5542.556058 


25.7772 


9.7 


-0.0461 


5542.748534 


25.6809 


11.3 


-0.0453 


5543.671498 


25.3436 


14.5 


-0.0548 


5557.628117 


25.5093 


10.1 


-0.0214 


5561.641322 


25.6004 


10.9 


-0.0256 


5562.627319 


25.4723 


9.6 


-0.0088 


5563.667964 


25.6645 


13.0 


0.0068 


5564.633499 


25.3992 


10.5 


-0.0214 


5565.645557 


25.6610 


11.8 


0.0061 



Table 1. CORALIE radial-velocity measurements for WASP-50 
(BS = bisector spans). 



tal abundances were determined from equivalent width measure- 
ments of several clean and unblended lines. A value for micro tur- 
bulence (^ t ) was determined from Fe i lines using the method of 
Magain (1984). The quoted error estimates include those given 
by the uncertainties in T e g, log g and £ t , as well as the scatter 
due to measurement and atomic data uncertainties. 

The projected stellar rotation velocity (vsin/*) was deter- 
mined by fitting the profiles of several unblended Fe i lines. A 
value for macroturbulence (v mac ) of 2.1 + 0.3 kms -1 was as- 
sumed, based on the tabulation by Gray (2008). An instrumental 
FWHM of 0. 1 1 + 0.01 A was determined from the telluric lines 
around 6300A. A best-fitting value of v sin i* = 2.6 ± 0.5 km s _1 
was obtained. 



Parameter 


Value 




5400 ± 100 K 




4.5 ±0.1 


ft 


0.8 + 0.2 km s -1 


v sin i, 


2.6 ±0.5 km s- 1 


[Fe/H] 


-0.12 ± 0.08 


[Na/H] 


0.00 ± 0.10 


[Mg/H] 


-0.06 ± 0.03 


[Si/H] 


-0.01 ± 0.09 


[Ca/H] 


-0.04 ± 0.08 


[Ti/H] 


0.00 ± 0.08 


[V/H] 


0.03 ±0.12 


[Cr/H] 


0.00 ± 0.08 


[Co/H] 


0.02 ± 0.04 


mm 


-0.07 ± 0.06 


log A(Li) 


< 0.6 ± 0.2 


Mass 


0.91 ± 0.07 M Q 


Radius 


0.88 ±0.11% 


Sp. Type 


G9 


Distance 


230 ± 40 pc 



Table 2. Stellar parameters of WASP-50 from spectroscopic 
analysis. 

Note: Mass and Radius estimate using the calibration of Torres 
et al. (2010). Spectral type estimated from T e $ using the table in 
Gray (2008). 



There is no significant detection of lithium in the spectra, 
with an equivalent width upper limit of 3mA, corresponding to 
an abundance upper limit of log A(Li) < 0.6 ± 0.2. The lack of 
lithium suggests an age ^ 0.6 Gyr, since similar levels of lithium 
depletion are found for stars of the same spectral types in the 
Hyades and Praesepe clusters (Sestito & Randich, 2005). 

The rotation rate (P rot = 16.85 ± 3.7 d) implied by the v sin i„ 
(assuming sin/* = 1) is in excellent agreement with the period 
measured from the WASP photometry, 16.3±0.5 days. This latter 
implies a young system, the deduced age being ~ 0.8 ± 0.4 Gy 
using the Barnes (2007) gyrochronology relation. 

The presence of Ca H+K emission (Fig. 6) indicates that 
WASP-50 is an active star, the deduced spectral index log R' HK 
being 4.6^ This value is similar to that found for the mag- 
netically active star WASP-41 (Maxted et al. 2011). Thus the 
spectroscopic activity indicates that the observed photometric 
variability of the star (Sect. 2.1) is due to the presence of spots. 



3 We note that a transformation for \ogR' HK measured with CORALIE 
spectra to a standard system does not yet exist. 
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Fig. 6. A region of the co-added CORAL IE spectrum of WASP-50 
showing emission in the core of the Call H and K lines. 



3.2. Global analysis 

We performed a global analysis of the follow-up photometry 
and the CORALIE RV measurements. The transit timing deduced 
from the analysis of the WASP photometry, 2454997.4996 ± 
0.0014 HJE0, was also used as input data to better constrain the 
orbital period. The analysis was performed using the adaptive 
MCMC algorithm described in detail by Gillon et al. (2009b; 
2010a; 2010b). To summarize, we simultaneously fitted the data, 
using for the photometry the transit model of Mandel & Agol 
(2002) multiplied by a baseline model consisting of a different 
second-order time polynomial for each of the four light curves, 
and a Keplerian orbital model for the RVs. 

The jump parameters^ in our MCMC analysis were: the 
planet/star area ratio (R p /R*) 2 , the transit width (from first to last 
contact) W, the parameter b' — a cos i p /R* (which is the transit 
impact parameter in case of a circular orbit), the orbital period P 
and time of minimum light To, the two parameters yfe cos a> and 
■\Je sin co where e is the orbital eccentricity and a> is the argument 
of periastron, and the parameter K2 = K Vl - e 2 P 1 ^, where K 
is the RV orbital semi-amplitude. We assumed a uniform prior 
distribution for all these jump parameters. The photometric base- 
line model parameters and the systemic radial velocity of the star 
were not actual jump parameters, they were determined by least- 
square minimization at each step of the MCMC. 

We assumed a quadratic limb-darkening law, and we allowed 
the quadratic coefficients u\ and «2 to float in our MCMC anal- 
ysis, using as jump parameters not these coefficients themselves 
but the combinations c\ — 2 x u\ + 112 and C2 = u\ - 2 x U2 to min- 
imize the correlation of the uncertainties (Holman et al. 2006). 
To obtain a limb-darkening solution consistent with theory, we 
decided to use normal prior distributions for u\ and U2 based on 
theoretical values and 1-cr errors interpolated in the tables by 
Claret (2000; 2004): u x = 0.29 ± 0.035 and u 2 = 0.29 ± 0.015 
for the / + z filter, and u x = 0.43 ± 0.025 and u 2 = 0.28 ± 0.015 
for the r filter. 

Our analysis was composed of five Markov chains of 10 5 
steps, the first 20% of each chain being considered as its burn-in 
phase and discarded. For each run the convergence of the five 

4 All the HJD presented in this work are based on the UTC time stan- 
dard. 

5 Jump parameters are the parameters that are randomly perturbed at 
each step of the MCMC. 



Markov chains was checked using the statistical test presented 
by Gelman and Rubin (1992). The correlated noise present in 
the light curves was taken into account as described by Gillon et 
al. (2009b). For the RVs a 'jitter' noise of 7.5 ms -1 was added 
quadratically to the error bars, to equalize the mean error with 
the rms of the best-fitting model residuals, 13 ms -1 . This jit- 
ter level is compatible with the activity level of the star (Wright 
2005). 

At each step of the Markov chains the stellar density deduced 
from the jump parameters, and values for and [Fe/H] drawn 
from the normal distributions deduced from our spectroscopic 
analysis (Sect. 3.1), were used as input for a modified version of 
the stellar-mass calibration law deduced by Torres et al. (2010) 
from well-constrained detached binary systems (see Appendix A 
for details). Using the resulting stellar mass the physical param- 
eters of the system were then deduced from the jump parameters 
at each MCMC step. 

Our measurement of the rotational period of the star P m 
made it possible to derive the posterior distribution for the stellar 
inclination /„. As in the analysis of Watson et al. (2010), at each 
step of the Markov chains we used the value of the stellar radius 
R„, and values for P rot and v sin /*drawn from the normal distri- 
butions JV(16.3,0.5 2 ) days and N(2.6, 0.5 2 ) kms -1 , to derive a 
value for sin z, from 



For sin i* > 1, sin i, was set to 1 and the inclination z* was set to 
90 deg. 

The parameters derived for the WASP-50 system can be 
found in Table 2. To assess the reliability of the deduced physical 
parameters we performed a stellar evolution modeling based on 
the code CLES (Scuflaire et al. 2008), using as input the stellar 
density deduced from our MCMC analysis and the effective tem- 
perature and metallicity deduced from our spectroscopic anal- 
ysis. The resulting stellar mass was 0.84 ± 0.05 M , in good 
agreement with the MCMC result, 0.89 + 0.08 M@. This stellar 
modeling led to an age of 7.0 ± 3.5 Gyr for the system. 

Fig. 7 shows the marginalized posterior distributions for 
the planet mass, radius, orbital inclination and eccentricity. As 
shown by the narrowness of these distributions, the physical pa- 
rameters of the planet are well constrained, thanks to the precise 
determination of the transit and orbital parameters made possi- 
ble by the high-SNR of our follow-up data (transit photometry 
and RVs). In particular, the orbit appears to be close to circular, 
any eccentricity greater than 0.05 being firmly discarded. 

4. Discussion 

WASP-50 b is a new planet slightly larger (1.15 ± 0.05 R Sap ) 
and more massive (1.47 ± 0.09 Mj up ) than Jupiter, in a close 
orbit (0.0295 + 0.0009 AU) around a moderately active G9V star 
(0.89 ± 0.08M o , 0.84 + O.O3fl ) at 230 ± 40 pc from the solar 
system. Taking into account the large irradiation received by the 
planet (8.7 + 1 x 10 8 erg s _1 cm -2 ), the measured planetary mass 
and size are in good agreement with basic models for irradiated 
planets (Fortney et al. 2007). Using these models our measured 
planet mass and radius suggest that the core of WASP-50b repre- 
sents only a few per cent, or at most 10 per cent of its total mass, 
depending on whether we take the age as ~ 1 Gy, as suggested 
by the rotation period of the star, or as older, as suggested by our 
stellar evolution modeling. Fig. 8 shows the location of WASP- 
50b in a mass-radius diagram, showing only planets with masses 
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Fig. 7. Marginalized posterior distributions for the planet mass, 
radius, orbital inclination and eccentricity resulting from our 
global MCMC analysis. 



between 1 and 2 Mj up . In this mass range the most extreme plan- 
ets are the very dense CoRoT-13b (Cabrera et al. 2010), and 
the highly bloated WASP-12b (Hebb et al. 2009; Chan et al. 
2011). WASP-50b appears to be a much more 'normal' planet 
with a density close to Jupiter's. We note however that the ex- 
tremely high density of CoRoT-13b was recently questioned by 
Southworth(2011). 

With a size larger than Jupiter's, a high irradiation, and a 
relatively small and infrared-bright (K=10) host star, WASP- 
50 b is a good target for near-IR thermal emission measurements 
with the occultation photometry technique, using Warm Spitzer 
at 3.6 and 4.5 yum (e.g. Deming et al. 2011) or ground-based 
instruments at shorter wavelengths (e.g. Gillon et al. 2009b). 
WASP-50 b orbits an active star (as does, for example, WASP- 
41 b, Maxted et al. 2010), making it able to test the proposition 
by Knutson et al. (2010) that the presence of thermal emission 
in the day-side atmosphere of a hot Jupiter is dependent of the 
chromospheric activity of its host star. This hypothesis suggests 
that an optical high-opacity source is present at high altitude in 
the dayside atmosphere of hot Jupiters and causes a thermal in- 
version, but that it is destroyed by the large UV flux produced by 
chromospherically active stars. This hypothesis is supported by 
emission measurements obtained so far, but it is still based on a 
small sample of systems. 

The values deduced in this work for log R' HK and the planet 
density p p are consistent with the correlation between these two 
parameters noticed by Hartman (2010). As mentioned above, 
stellar evolution models and the measured density and spectro- 
scopic parameters of WASP-50 point to an age of 7 ± 3.5 Gy, 
translating into an expected rotational period of 50 ± 15 days un- 
der the empirical relationship of Barnes (2007), calibrated on a 
large sample of cluster and field stars. Our measured rotation pe- 
riod (16.3 ± 0.5 days) is in poor agreement with this expectation 
and suggests a younger age. In the context of a \ogR' HK -\ogg p 
correlation, this discrepancy could be explained by the action of 
the planet itself on the star. In this scheme, hot Jupiters orbit- 
ing moderately cool stars like WASP-50 would tidally spin-up 
the stellar convective envelope and tilt its rotation axis, the stel- 
lar magnetic braking preventing the system from reaching spin- 
orbit synchronization, in agreement with the observations that 
stars hosting a hot Jupiter tend to rotate faster (Pont 2009), and 
that hot Jupiters with a high stellar obliquity are preferentially 
found around hot stars having no convective envelope (Winn et 
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Density p, [p Q ] 
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4.537 ± 0.022 


Mass M, [M Q ] 


892+ 080 


Radius R, [R e ] 


0.843 ±0.031 


sin i. 
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Deduced planet parameters 


RV«[m s" 1 ] 


256.6 ± 4.4 


R P /R, 


0.1404 ±0.0013 


b tr [R,] 


0.687t™ 


b oc [R,] 


0.691!°™ 


T oc - 2450000 [HJD] 


5561.549 ±0.010 


Orbital semi-major axis a [AU] 


0.02945 ± 0.00085 


a/R, 


751+ 0.17 


Orbital inclination i p [deg] 


84.74 ± 0.24 


Orbital eccentricity e 


0.009!™" 


Argument of periastron a> [deg] 


44+62 


Equilibrium temperature T eq [K]" 


1393 ± 30 


Density p p [p Jup ] 


958+ 095 


Surface gravity log g p [cgs] 


3.439 ± 0.025 


Mass M p [M Jup ] 


1.468_ 0086 


Radius R p [R lup ] 


1.153 ±0.048 



Table 3. Median and l-cr limits of the posterior distributions ob- 
tained for the WASP-50 system derived from our MCMC analy- 
sis. "Assuming A=0 and F=\ . 



al. 2010b). Moreover, we note that our deduced stellar inclina- 
tion i, (84^^ deg), although imprecise due to uncertainty in the 
v sin determination, argues against a large obliquity. 
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et al. 2010, CIO) and the ones recently reported by South worth 
(2011, Sll). 
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Appendix A: Stellar calibration law 

Using a large sample of well constrained detached binaries, 
Torres et al. (2010; hereafter T10) showed that accurate stellar 
masses and radii could be deduced from an empirical law relat- 
ing these two parameters to the measured spectroscopic param- 
eters T^, logg, and [Fe/H]. Enoch et al. (2010) showed that 
this law could be advantageously used in the analysis of tran- 
siting planets data with the MCMC algorithm, provided that the 
stellar surface gravity was replaced by the stellar density in the 
calibration law. Indeed, the stellar density is directly and pre- 
cisely constrained by high-quality transit photometry, making it 
possible to obtain better stellar parameters than by using the stel- 
lar surface gravity deduced from a spectroscopic analysis (Winn 
2010a). Under this condition, the calibration law makes it possi- 
ble to deduce the mass of the star, and thus the physical parame- 
ters of the system, at each step of the MCMC, without the need 
for a concomitant or subsequent stellar evolution modeling. The 
main limitation of this approach is that the T10 law is poorly 
constrained for stars below 0.6 M and for pre-main-sequence 
stars. 

The calibration law of Enoch et al. (2010) relating M* to p*, 
T e ff and [Fe/H] is 



log M, = ci\ + aqX + CI3X 2 + a4(logp«) + 

a 5 (logp*) 2 + a 6 (logp») 3 + a 7 [Fe/H], 



(A.l) 



where X is log T e g -4.1. 

We performed a linear regression analysis with a Singular 
Value Decomposition algorithm (Press et al. 1992) to derive the 
best-fit values for the a\.j coefficients, using as data the tabulated 
values presented by T10 for 19 binary systems for which the 
metallicity is known. We added another very well characterized 
star to these data: the Sun, assuming conservatively a relative 
error of 10~ 3 on its mass and radius, 0.01 dex on its metallicity 
and 10 K on its effective temperature (5777 K, Cox 2000). 1-cr 
errors were derived through a Monte Carlo analysis. The fitted 
and error values for the a\j coefficients are shown in Table A.l. 
The resulting scatter in fitted-minus-measured values is 0.084 
M (0.023 in logM,). 
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Table A.l. Coefficients for stellar mass calibration law. 





0.450 ± 0.013 


a 2 


1.421 ± 0.057 


03 


0.33 ±0.19 




-0.082 ±0.016 


a 5 


0.032 ±0.017 


a 6 


0.0026 ± 0.0037 


a-, 


0.075 ± 0.012 



To take into account the internal uncertainties of the calibra- 
tion law, the a i;7 coefficients of Eq. A. 1 are not fixed during our 
MCMC analysis but are drawn randomly from the normal distri- 
butions shown in Table A.l at each step of the Markov chains. 



